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I. INTRODUCTION 

Single deformable objects such as droplets of one liquid immersed in another liquid 
fluctuate in shape in response to random external stirring pL]. The purpose of this paper is 



to investigate the effect of interaction between such deformable objects [2] on the way they 
deform due to fluctuations of the velocity field in the fluid in which they are immersed. 
Numerous authors have studied the fluctuations and diffusion of single deformable objects 
due .„ the..al agitation g fl fl fl S B B Q. Neve.h*.. tWe a. c,ea.,y otW 
ways in which systems are agitated. In industrial and biological environments, the host 
liquid is often stirred, shaken or pumped in ways which are very different from thermal 
agitation. The list of examples is not restricted to artificial processes. It also includes 
natural processes such as Brownian motion of small beads induced by the collective motion 
of bacteria [11] and nano-scale mechanical fluctuations of the red blood cell surface that 
have been measured and shown to depend strongly on the biochemical environment and 



not only on temperature |l2l . Il3l . Il4l . Il5l | . For this reason the external velocity field agitating 
the system is taken to be more general than that corresponding just to thermal motion. 
The article provides thus the general equations describing the effect of a finite density of 
deformable objects on the diffusion and shape fluctuations of a single object, to linear order 
in the density. 



The system considered has the following properties, 
(a) The deformable objects are fluid, in the sense that the velocity field is well defined 
everywhere (both inside and outside the object). No slip and no penetration conditions are 
assumed at the interface of the deformable object. Hence, each surface element moves with 
the velocity of the flow at its position. In addition, both the objects and the host fluid are 
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incompressible. The objects are characterized by an energy that depends on their shape (i.e. 
changing the orientation or switching places of two surface particles while keeping the shape 
constant does not change the energy). The main example considered in this paper is surface 
tension Ig], . However the description can be extended to Helfrich bending energy [3] 
and other cases where the shape of minimum energy is nearly spherical. Deformation of the 
shape changes the energy, exerts a force density on the liquid and therefore generates an 
additional velocity field, denoted by v^. 

(b) The hydrodynamic equations of the host liquid are linear in the velocity ( i.e. a velocity 
field induced by several sources is equal to the sum of the velocity fields induced by each 
source separately). For instance, if the flow is governed by the Navier Stokes equation, then 
the linearity implies that the Reynolds number is small and that the Stokes approximation 
is applicable. The actual velocity field is the sum of the imposed velocity field , Vext (the 
velocity field that would have existed if the objects were absent), the velocity field induced 
by the deformations of the object under consideration, v^, and Vr which is the velocity field 
created by the rest of the deformable objects, 

V =Vext + V^+Vr. (l) 

(c) The external velocity field is assumed to be random. The correlations are assumed 
depend only on distance and time difference. Furthermore, the dependence on the time 
difference is taken to be extremely short ranged (Dirac 6 function in the time difference). In 
principle, equations for the dependence of the shape correlations on the density of deformable 
objects can be worked out for any dependence of the velocity correlations on time. These are 
very complicated, however, and the above choice of the dependence of the external velocity 
correlations on time simplifies matters considerably and is certainly realistic in many cases. 
It is important to note that the results obtained here are not used to determine the external 
velocity correlations. Those correlations are just taken as a given input. For example, in the 
special case of thermal agitation the velocity correlation was calculated from first principles 



10|] and only then used to calculate the diffusion constant and deformation characteristics 
of a deformable object immersed in the liquid. In addition the external velocity is assumed 
to be small enough to allow the body to remain almost spherical. 

Since we assume small deviations from the spherical shape it is only natural to describe 
the surface shape of the objects using spherical harmonics. Consider a spherical body which 
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is moving and is slightly deformed. The equation 

^ + f{Q,t) -1 = (2) 

defines its surface, yielding for each spatial direction, Q, the distance, p = |r — ro|, of 
the surface from the centre of the body, tq. R is the radius of the undeformed sphere. 
The deformation function, f{Q,t), defines the shape and may be expanded in spherical 
harmonics, f{^,t) = J2'iliYlm=-i fi-m{'t)Yim{^) (clearly the loo term can be absorbed 
in the definition of R). The goal is to obtain the correlations between the deformation 
coefficients, fim{t)- The centre of the object, ro, is chosen to be the point around which the 
deformation coefficients with / = 1 vanish: /i^ = 0. A different definition of the centre will 
introduce three different equations for the deformation coefficients with 1 = 1. These are 
not interesting, as far as the shape is concerned, since in the first order of the deformation 



the spherical harmonics with / = 1 describe a rigid translation of the object 20|, l2l[ |. 



The random velocity field and the effect of the interaction between objects induce fluctu- 
ations in the values of the deformation coefficients describing any of the objects. Consider 
the autocorrelation of the deformation coefficient of the i'th object {flmi^)fl-mi~'^))y 
where flmiy^) and flmil^) represents the Fourier transforms (FT) of the deformation coeffi- 
cients of the i'th object in respect to time. The autocorrelation is expanded in orders of n. 
To first order in n it is given by 

{flm{^)fl-m{-i^)) = ^/,0(^,^) + Ql,l{uj,u) . n (3) 

The first term on the right hand side of the equation (3) above gives the shape correlations 
of a single object , that has been described previously [l| . The second term represents the 
correction to the shape correlations due to a small but finite density of deformable objects. 
The purpose of this article is to obtain that correction. 

The paper is organized as follows. Section [TTl deals with a single deformable object in 
random fiow. The aim of this section is to introduce the basic definitions of fiow and present 
the zero-order terms in the expansion of the shape correlations in the density of objects n. 
In section IIIII the first order terms of the shape correlation functions are derived. Section 
IIVI deals with the special case of identical droplets that are governed by surface tension in a 
random fiow. In order to improve readability. Parts of the derivation of the first order terms 



and the velocity field induced by a deformable body governed by surface tension are left to 
the appendix. 



II. A SINGLE OBJECT IN RANDOM FLOW 

The response of a single object immersed in a host liquid to an external random flow has 
been described in previous work fl]. The results are sketchily repeated here for the benefit 
of the reader as the general equations obtained here are to be exploited in the next section 
by replacing the external velocity field by the velocity field seen by the object when a finite 
density of deformable objects is immersed in the liquid. The latter velocity field is the sum of 
the imposed velocity field and the velocity fields induced by all the other deformable objects. 
The correlations of the deformation coefficients as well as diffusion constant of the centre 
will be obtained here in terms of the correlations of the external velocity. The diffusion 
constant will be needed in the next section to obtain the first order correction in the density 
of the shape correlations. 

The no-slip and no-penetration conditions yield an equation of evolution for the deformation 
coefficients {2^, 

The effect of the velocity field induced by the deformable object itself is represented by 
the second term on the left hand side above. The A;'s characterize the way in which a 
deformation with definite / decays to zero in the absence of an external velocity and other 
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objects. Different physical systems are characterized by different sets of A; 
The term on the right hand side, Qim, is given by 

Qi-- = ^j ■ [^e.t - ri,] Y*jn)] , (5) 

where the external velocity field, Vext, is evaluated on the undeformed body and is a unit 
vector in the direction of the spatial angle Q (for further detail see [l^). The velocity of the 
centre vq contributes only to Qi^m- In addition, the definition of the centre implies that 

Ql,m = 0, (6) 

for all m. 

It is convenient to write the correlation function of the external velocity field in momentum 



space. This is so because the random velocity field is transversal when the fiuid is incom- 
pressible. Consequently, in real space, the fiow must always be correlated in a very complex 
way. On the other hand, in momentum space, the transversal part of a general field is easily 
obtained: 

QiQj 



(7) 



where u is a general vector field and i and j denote Cartesian components. The bracketed 
term is the projection operator that removes the longitudinal part of u, and therefore yields 
a general transverse velocity field Vext- Next, the correlations of the external velocity are 
easily expressed using the correlations of the general field u. 

{u'{q,t,)u"'{p,t2)) = 6i^6{q + p)(j){q,\t2 - (8) 

where is a general function of q and the time difference (with the only limitation that 
its Fourier transform in the time difference is non- negative). As was mentioned before, this 
investigation is restricted to cases where the external velocity is uncorrelated in time, 

0(g,t)=0(g)5(t), (9) 

and where the mean of the velocity field vanishes, (^u\q,t)'^ = 0. 

Jsing the above, the diffusion coefficient of the centre of the deformable object is obtained 



2l|, 



D 



dq 4>iq)J\qR), 



(10) 



where J'{x) = jo{x) + j2{x) and jn{x) is the spherical Bessel function of order n. 
The correlations of the deformation coefficients, fi^rn, |l| are given by 



{flm{t)fl'm'{t + At))^_ 



Xi\At\ 



2X, 



I' I '-'m'—rm 



(11) 



where 



Qw I dn I dn' I d?q Y*in)Y;,in') 



(Mi 

g2 



(12) 



The contribution to the autocorrelation is presented here for completeness after being Fourier 
transformed in time, 

6{uj — Lj) 



47r2(A2 + tu2) 



(13) 
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III. FIRST ORDER CORRECTION 



The aim of this section is to obtain the correction to the autocorrelation of the object's 
shape, which is hnear in the density of objects n. The full autocorrelation of deformation 
coefficients is obtained by replacing Vext in equation (jlj) by the sum of the external velocity 
and the velocity induced by the other objects, 

G\q2, -uj - uj-i) ^i) + ^i)) (wL(g2, ^2) + (g2, ^2)) 

In this and in all the following equations, the Einstein summation convention is applied to 
the Cartesian components a,/?. In addition, G^{q,uj) appearing in the equation above, is 
defined as the temporal FT of 

G*(g,t)=exp(^g-r'oW)- (15) 

The first order correction in the density has two contributions. The first contribution is 
obtained by neglecting Vr in both brackets on the right hand side of eq. f|T4l) but taking the 
G"s to first order in n. This results in a contribution given by 

ng\l}{u;,u) = ~n5{uj-^)^j^{2^)-^R-^ j dn, j dn^ j d'qY*^{n^)Y*_^in2) (16) 

q^ J 6 



In this general form JF = JFq + jFi ■ n is the total mean square displacement (MSD) of the 
center of mass. The above contribution vanishes, however, for cases where the bare velocity 
field is uncorrelated in time, eq. (Q. Consider the integral over t on the right hand side of 
eq. ([9]). The Dirac delta function in (f){q,t), eq. ([9]), sets t = 0. Since for t = the MSD 
must vanish in any order of the expansion in n, JFi(O) = and the right hand side of eq. 
f ll6p vanishes. 

The second contribution arises from those terms in (fT^ linear in Vr, taking the MSD to 
zero order in n. In the linear approximation and for small deformations, the velocity field 
created by the deformable objects can be written as a sum of prefactors, Oj^ that multiply 
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the deformation modes of the j'th object. Considering a set of identical objects, it is obvious 
that Oj^ is identical for all j. Thus, the superscript is dropped and tv is written as 



jy^i l,m 



(17) 



The prefactors O/ ^ are model dependent and are calculated in appendix [B] for deformable 
objects that are governed by surface tension. Using the above expression for Vr, the second 
contribution is given by 



(cj,cD) = -27^<( 6{u-u) 



1 



1 



(18) 



I'm' 
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g| iuj + \ii + qfD 



0(^2 



where TZ{x} denote the real part of x (for detailed derivation of equation (fT8|) see appendix 



Once Oi m and 



2d . I22I . l23l | of a specific system and the correlations of the bare 



velocity field are known, the correction to the shape correlations can be calculated using the 
above equations. In what follows, the actual use of the general equations is demonstrated 
for a specific system of deformable objects that are governed by surface tension. 



IV. DROPLETS GOVERNED BY SURFACE TENSION 

To determine the shape correlations we need 0(g) that describes the effect of external 
agents on the system. We need the set of the A;'s describing the decay of a deformation of 
angular momentum / of a given membrane in the absence of the bare velocity field. We also 
need the 0/ ^'s that describe how the velocity field in the liquid responds to the deformation 
of a single object. The above quantities depend on the properties of the deformable objects 
and of the liquid in which they are immersed. 

The system to be considered in the following is a system of membranes governed by surface 
tension. Namely, the energy of the membrane Us is given by Us = XS where 5* is the surface 
area. The viscosity, t], is assumed to be uniform inside and outside the objects. (The 
qualitative behaviour for different viscosities is not changed for a finite reasonable range of 
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viscosity ratios 2J] and using identical viscosities eliminates the boundary conditions, thus 



simplifying the calculations considerably). For that case the A^'s were derived in the past 

_ A (/ + 2)(/ + l)/(/-l) 

The 0/,m's, for the same system, are needed only for large distances from the centre of the 
object inducing the velocity by its deformations, because the density of deformable objects 
is low. The leading non trivial contribution to the Oi^m is calculated in appendix [B] and is 
given by 

Ol,m{f^ = 6l,2^ (^^^ ^Y.^mim, (20) 

where Cl is the unit vector in the direction of f. 

Now, the integration over Qi, Q2 and ^3 in eq. (ITSl) can be easily done using the partial 
wave expansion, 

00 I 

^-rm = J2J2 {-^y^^JlmyLwyln.m■ (21) 

/=0 m=-l 

These integrations produce a long but finite set of terms that is not presented here because 
of its length and complexity. To facilitate the integration over the g's, the following tactics 
is used. The above expression is split into three parts by defining the following Fourier 
transforms: 



d^i / d'q, . Ol^mYU^iWi-—^. (22) 

J iuj + X2 + qtD (27r)2 



B^Mr2) = j dn^j dQs j d%e''^^<^^+^'^^{6,,^ - ^)0(92) (23) 



(^2 ) ^2*m' (^3) ^2 ^3 ' 



(2vr) 



2 



and 



d'qJs,^^-l) . (24) 

^ ^ [ZTTj 2 
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Combining the above together, the first correction to the shape correlation for a dilute 
system is obtained by changing the order of integration and performing first the integration 
over the gr's, 



Note that since the system is supposed to be invariant under rotations, the structure 
factor can depend only on the absolute value of q and C must depend only on the absolute 
value of r. 

An additional property of the shape correlations is that the expression: 
2iTT {fi,mi^)flm*i^'^)) transforms as a scalar. Moreover, due to the rotational symme- 
try, the shape correlation does not depend on m and therefore, 

{flMfUi-^)) = ^ E {flm'i^)fUi-^)) ■ (26) 

m' 

In this way the coordinate system is easily rotated without loss of generality. 

Last, note that the expression for C{r) is, up to a prefactor of n(27r)~i, no other than 
the pair distribution function. The pair distribution function C is chosen, as a good approx- 
imation in the dilute regime, to be the pair distribution function of hard spheres system in 
the dilute limit, 

^ , f if r < 2i? 
C(r) = { (27) 
l^n^ if r > 2R. 

A straightforward but tedious derivation, shows that the shape correlations for all I except 
I — 2 vanish and for Z = 2 it is given by the following expression: 



(3cos(qR)qR - 3 sm( qR) +sm(qR)q'^ R"^) , ^ . , ^ . , ^ 

dq2(p{q2) sm{qr) sm{q2r) 



(lu + X2 + q^D)q^ jQ 
(-9 cos{q2R)q2R + 9 sm{q2R) - 4 sm{q2R)q^R^ + cos{q2R)qiR^f 

The integration over q can be done analytically but will not be presented here due to the 
length of the expression. Note also that the integration over r must be done last. Once is 




FIG. 1: The correlations of the deformation coefficient with / = 2, {f2mi^)f2mi~'^)*)^ t° zeroth 
order (dotted Une) and first order (continuous hne) in the density of objects. The units along the 
y axis are relative and the units along the x axis uj are normalized by the decay rate A2. 

known, the result can be calculated analytically for special cases or numerically for others. 
A specific example for the use of eq. ( !28l) is depicted in fig. [T], for a correlation function that 
has the form of the Yukawa potential. The Fourier transform of the correlation function is 
given by 

4>{<l2)=<Po/{a' + ql), (29) 

where a = is used in this example. In addition ^^^^ = 1 is used. The correlation of 
the deformation coefficients, //^, with / = 2 is depicted to zeroth order (dotted line) and 
first order (continuous line) in the density of objects. The correlation of the deformation 
coefficients must have the form, 

where F is a general function that depends on I and u. In general, the correction that is 
linear in density, nQi, must be small relatively to the zeroth order term Qq. In this example 
however, the density n is chosen to be large enough to observe changes. The correlations 
of deformation coefficients with Z 7^ 2 do not change to first order in the density of the 
deformable objects. As can be seen, the deformation modes with / = 2 are suppressed by 
the interaction between the droplets at low frequencies, u < X2. This is expected due to the 
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retarded response of each droplet to the external velocity field. This decay rate A2 introduces 
a new time scale that controls the decay of fiuctuations produced by the external field. At 
low frequencies, lower than the time it takes the deformation modes to decay; the velocity 
field, induced by neighbouring droplets, responds effectively to the external velocity field 
and thus decreases the fi^m correlation at low frequencies. 

APPENDIX A: ON THE CALCULATION OF THE DEFORMATION CORRE- 
LATIONS 

This appendix derives the contribution of the terms involving t>"^j(g, c<Ji)f^(g2, ^^2) and 

t>"(g, ciJi)t>f^^(g2, 1-^2) in eq. f|T4|) (i.e. Q^'^i). These two terms are complex conjugates of each 

(2) 

other. Thus only the first term will be considered and the correction Q^^ is given by twice 
the real part of the answer. In order to keep the expressions to first order in the deformation 
and density of objects, expressions must be kept linear in Vr and fi^m- First, the average is 
broken into two parts, the velocity correlation and the average over expressions containing 



fo [1^ ( This approximation was justified and used a number of times in the past 



Id, 



25|) 



j^i l',m' 



The Use of the expression for the correlation of the bare velocity ,eq. ([8]), and the 
definition of G\ eq (fT5l) . yield 
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f d^, /<i^2V;:„(Ji)y,:_„(SJ2)e*'»''e««.«eS;«lEEO?.".'(9) 

j^i l',m' 



G\q,uj - uJi)G\q2, -Cj - uj2)G\-q,uJi - uJz)G\-q2,0J^ + UJ2) 



where 



G\q,u-uJi)G\q2,-uj-uj2)G^{-q,uJi-uJz)&{-q2,uj^ + uJ2)j = (A3) 
exp(zg • (rl,(ti) - rSlts))) exp(^g2 ■ (r-*o(t2) - r^oiU))))- 



(27r) 



The integration over loi can be done using J dcoie'^^'^^*'^ — 27rS{ti — ts). Decouphng the 
average into a product of averages one finds that 



J2 ( exp (iq- (r- (ti) - f^ih))) exp (zgl • {flit2) - f^U))) ) (A4) 

= ^ (exp (zq. {flih) - fl{t2))) ) (exp ( - zq- {rUh) - f^iU))) ' 

(exp (t{q + q2) ■ (^0(^2) - fl{U))) ) (exp ( - t{q + gl) ■ {rj^iU) - ^0(^4))) ^ 
= exp ( - ^T{tr - t2) - - U) - ^^^^Ht2 - U)) 

^ ^ (exp (i{q + q2) ■ fl{U)) exp ( - i(g + ^2) • r^(i4)) ) , 

i jjti 

where N is the number of deformable objects in the system. The reason for the decouphng 
of exponents is that the driving external velocity correlation decays both in time and with 
distance. Because of the low density and the short range repulsion (as of hard spheres in 
our approximate description to the first order in the density of objects) the exponents can 
be only weakly correlated. 

The addition and subtraction of the i'th object to the sum over j results in 



e 6 



J2 { exp(i? • (rl(ii) - fiim eMiQ2 ■ (rt(t2) - ^0(^4)))) = (A5) 
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By combining all the above together, using variable transformation and performing the 
integrations over tua, uj2 and t4, the shape correlation are derived. 



^gV,^) = -27^|5(^ - jdnj dn, J d\ j d% (A6) 



l',m' 



7 

^3 



7 /3 



oo /"OO 



e-^-<'=>(s,,„-l)}, 



from which (fTSl) is obtained. 



APPENDIX B: THE VELOCITY FIELD GENERATED BY DEFORMATION 
FOR THE SURFACE TENSION CASE 

Consider a liquid droplet governed by surface tension that is immersed in a host liquid. 
Assume that the viscosities inside the droplet and in the host liquid are equal. The defor- 
mation of the shape of the object changes its energy and in response induces a force density 
that acts on the fluid. The force density creates in its turn an additional velocity field, 
denoted here as v^, 

y^(f) = i y S{f- ■ F{r')dh\ (Bl) 
where F is the force density created by the object and S is the Oseen tensor that is given 



by 



= + (B2) 



It is easy to calculate the force density using simple field theory. Let %l){r) be a three 
dimensional scalar field, defined everywhere in such a way that the equation %l){r) = 
describes the surface of the object {2!, U^- The gradient of ip is assumed to exist and not to 
vanish in the vicinity of ipl-r) =0. Under the additional assupmtion that the deformation of 
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the object do not produce over-hangs,-^ (r^ is written using the deformation function, f{fl), 
given in eq. ([2]). 



(B3) 



In 



20| the force density created by a deformed objects governed by surface tension is 



given by. 



F(f) = -A(V ■ n)6{ij{r))Vij{f) 



(B4) 



where n = is a unit vector in the direction normal to the surface of the deformable 
object. 

In addition, the following symbols are used for the angular parts of the gradient and 
Laplacian, 



V = 0— + (f- 



and 



36 ^?,ui{e)d^ 



sm^e')— 1 + 



(B5) 



(B6) 



sin(^)9^ V""^'"'^^y ' sin2(e)9(/?2' 

where 9 and are unit vectors of 9 and (p. The velocity field is obtained as a function of 
/(fi', t) by the use of equations fIBip and fIBSp . 



A 



1 ^' + T^ 



- / S{r~R{l~f)n') 

7] 



R 



2 , 

1-/ ^ 



1 + 



(1-/)= 



l-f 



2 JO/ 



;i-/)w. 



2 fi + il 

U + (l_/)2. 

Keeping the above expression to the first order of the deformation / results with 



(B7) 



XR 



dVt' 



§ + 2^X-^^RX + 6^1p^Rx)m') (B8) 



2Vf{n') , 2X-{Vfm) 



X 

Q' X -Q' 

X ^ X3 



X3 



X 
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— * A 

where X = r — RVt' . f is replaced by its expansion in spherical harmonics f{Q,t) = 
Z]|^2 ^m=-i /'.™(^)^'.™(^)- "^^^ spherical harmonics are eigenvalues of V^. In contrast, V 
mixes different harmonics. The density of objects is assumed to be low and therefore the 
distance r to the closest droplet is typically much larger than the radius of the droplet R. 
Hence, the velocity is expanded to the first nontrivial order in ^. In addition, a rotated 
coordinate system is used in which r = rz', where z' is the z direction in the rotated 
coordinate system. In that system, 

A /rV 2 



v4rz',t) = -[-] -j=f2,o^', (B9) 
V \ r J VStt ' 



where /20 is the deformation coefficient of 1^20 rotated coordinate system. The 

velocity field in a general direction is given by a rotation of the coordinate system. The 
transformation under rotation of the spherical harmonics (and thus of /2,o) is given by the 
addition theorem, 

47r ^ 

^n(cos(7)) = Yl Y,,,U0uVi)YZm{02,^2), (BIO) 

m=—n 

where 7 is the angle between the directions Cli and Q2 that corresponds to {6i,{pi) and 
{02, '^2)- Note that the calculation of f^Q involves integration over all angles that produces 
a dependence of the induced velocity on the angles at the original coordinate system. Thus 
the induced velocity in a general direction is given by. 



^ ^ E y2,rnmf2,mm (Bll) 

' ^ ^ m=-2 

where f is a unit vector in the direction of r. Oi^m is obtained by comparing eq. (fT7|) with 
eq. mH), 



Oi,m = Si,2-[-] -Y2,U^)f. (B12) 
T] \ r J 5 

Hence, the only spherical harmonics modes that contribute to the velocity field, far from 
the object, are the terms with / = 2. 
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